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Abstract 

Abstract: Bone surface mapping technique is proposed on the bases of two kinds of uniqueness of bone 
in vivo, (i) magnitude of the principal moments of inertia, (ii) the direction cosines of principal axes of 
inertia relative to inertia reference frame. We choose the principal axes of inertia as the bone coordinate 
system axes. The geographical marks such as the prime meridian of the bone in vivo are defined and 
methods such as tomographic reconstruction and boundary development are employed so that the surface 
of bone in vivo can be mapped. Experimental results show that the surface mapping technique can both 
reflect the shape and help study the surface changes of bone in vivo. The prospect of such research into 
the surface shape and changing laws of organ, tissue or cell will be promising. 

Introduction 

The shape of bone is the adaptive result of bone in the mechanical and physiological environment [IH4| . 
A map, the visual representation of our real world symbol model, can reveal not only the spatial structure 
properties of an object but also the changes in time series [5HZ]. The mapping of bone surface, therefore, 
is used as an approach to study the adaptability of bone morphology. Mapping technique has become 
an even more powerful and useful method to do scientific research. Mapping and flattening techniques 
have also been widely used in medical research [81-113] . They are both concerned with the development 
methods. The flattening technique develops the three-dimensional object to a two-dimensional one [14] 
while the mapping technique plays an essential role in interpreting the surface structure of an object |15j . 
The advancements and improvements of three-dimensional imaging of bone in vivo [16H18] have brought 
better data collection methods of bone surface, but the mapping techniques have not been systematically 
explored with satisfactory results. 

To map the bone surface, some geographic marks and geographic coordinates such as the bone surface 
prime meridian, the equator line or the contour should be identified pTS] . They should be defined on the 
bone's coordinate system. In order to study the bone's changes caused by the external factors, it is 
fundamental to set up a bone coordinate system when mapping the bone surface. In our study, we 
have set up a coordinate system with uniqueness on the principal axes of inertia of the bone in vivo. 
The coordinate system is thus employed to determine the prime meridian, and the average radius of the 
tomographic boundary is employed to determine the contour to map the bone surface. This means that 
the bone surface mapping technique could present an alternate approach to study the bone's morphology. 



Standardized coordinate system of bone 

The CT image of bone in vivo can generate the principal moments of inertia and the direction cosines |20) . 
Suppose the moment of inertia of the bone's tissue relative to its center of mass is a constant; then 
the magnitude of the bone's principal moments of inertia will be determined by its shape and mass 
distribution. The inertia tensor suggests that we can always find a group of coordinate systems where 
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three products of inertia will be nil at the same time [2TJ[22]. The magnitudes of these three principal 
moments of inertia could come out with three results: (i) all three are equal; (ii) two out of three are 
equal and (iii) each one is different from the other. When the object is homogeneous, in the first case, it 
is a sphere; in the second, an ellipsoid, a cube, a cylinder or a rectangular. When the bone's location and 
orientation relative to inertia reference frame (i.e. the coordinate system is defined by the CT coordinate 
system) are fixed, in the first case, there are numerous principal axes of inertia, while in the second 
case, the orientation of one principal axis can be determined, but not the other two. Therefore, in the 
first and second case, the principal moments of inertia have nothing to do with their principal axes of 
inertia. In the third case, however, the direction cosines of the bone's principal axis of inertia relative to 
inertia reference frame are unique, which means a one-to-one corresponding relation between the direction 
cosines of the bone's principal axes of inertia and its shape. 

If the bone is defined as a collection of elements of volume AV, the element's position can be repre- 
sented by (xgi, Doi, z i) (where o is located in the center of mass ), then the magnitude of the principal 
moments of inertia can be represented respectively by: 

4 = E(^ + 4) PiAv 

Iy = E (xli + 4) PiAV (1) 

i* = Y, + vl) pi&v, 

where p is the density, AV = AxAyAz, Ax and Ay are the pixel sizes and Az the layer distances of CT 
images. 

Let the angular displacements of the bone that take turns to rotate around axes x, y, z be a, /3, 7. 
According to Equations [TJ we can set up the following equation: 

(Iy - I z ) a = Y^[(VoiCOsa - z ol sina) 2 - (y oi sma + z oi cosa) 2 ]piAV 

(I x ~Iz)p = J2i( x ot cos p + z oi sin /3) 2 - (x ol sin/3™ z oi cos (3) 2 }piAV (2) 

(4 -Iy) 7 = J2i( x oi COS7 - y oi sin7) 2 - (x oi sin 7 + y oi cos7) 2 ]p i AV 

Next, differentiate Eq. [2] and let 

d(ly-I s ) a d{l X ~I Z ) p d(l X -Iy) 7 

da dp g?7 

which will generate the following equation: 

P = - arCtan U(^g(S.A,) ) (3) 

T^^4^)lt^) {c) - 

It is apparent that when and ONLY when I x ^ I y I z will Eq. [3]have a set of solutions. J] UiZiPi AV, 
Y^, ZiZiPiAV and ^ xiyipiAV are three products of inertia of inertia tensor. Within the range of [0, it], 
according to Eqs. [2]and[3j the limited rotations can always turn three products of inertia into zero at the 
same time. 

The bone's shape is asymmetrical and its structure is anisotropic [23H27] . Eq. [3] exposes the direction 
cosines of principal axes of inertia relative to the inertia reference frame of the bone characterization are 
unique. As a result, the coordinate system set upon the principal moment of inertia can not only depict 
the position and orientation of the bone, but also verify the bone surface shape and its changes when 
making a quantitative analysis. Eqs. [2] and [3] also suggest that we can set up a coordinate system whose 
coordinate origin is located arbitrarily at the center of mass of the bone. After limited rotations, the 
coordinate axes is positioned on the principal axes of inertia. 
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Mapping of the Bone surface 



CT scanning simplifies the bone as a collection of elements of volume AV (different densities). When 
performing an isotropic scanning (where pixel size must be the same as the layer distance), the volume 
of AV is a constant. The position of AV relative to the center of mass differs from one another. When 
the result of Eq. [3] is replaced for that in Eq. [2] accordingly and when the bone's coordinate system is 
positioned on the principal axes of inertia, the rotation will change the original CT image. It is necessary, 
then, to reconstruct the new image, which can be performed by the following equation: 



.r 



trunc 
Hi = trunc ' J 
Zj = trunc 



x i—m'm(x i) 
Ax 



Ay 
-min(z i ) 



Az 



Ax + min(x i) 
Ay + min(y oi ) 
Az + min(z i), 



(4) 



where (x i,y i, z Q i) stands for the position of AV after rotation, (xi,yi,Zi) for that of the reconstructed 
tomogram and truncQ for a function that truncates a number to an integer by removing the fractional 
part of the number. To keep the CT images isotropic, Ad x = Ad y = Ad z is defined in Eq. 01 and its 
pixal size and layer distance are kept the same of those of the original image. According to Eq. 21 when 
keeping the AV to be in a cube, the new CT image after rotation remains to be closed and continuous. 

CT scanning divides the bone surface into a collection of the tomographic image boundaries. In this 
way, the mapping of the bone surface has become an issue to develop the tomographic boundary, making 
it necessary to detect and draw the tomographic boundary. Before a new CT scanning, the equipment is 
reset, i.e. the gray value of the air is set as zero. The scanned tomographic images of bone are processed 
by Eq. 2J and their boundaries are drawn by the following equation: 



p(x,y)z 



p(x,y) 2 
p(x,y)z 

other, 



p{x,y)z 
p(x,y)z 
p(x,y)z 
p(x,y)z 



> 0,p(x+l,y) z 

> 0,p(x - l,y) 2 

> 0,p(x,y+l) z 

> 0,p(x,y - 1) 2 



(5) 



where z is the number of sequence of tomogram, (x,y) z the position of AV relative to the tomographic 
center of mass and p(x,y) z the density of (x, y) z . 

When mapping the bone surface, the bone will be "cut" - from a cylindrical surface to a rectangle, 
or a rhombus. However it is cut, its ultimate area would be the same. But when it is a rectangle, there 
is only one. How to cut the bone into a rectangle? The two principal axes of inertia (the minimal and 
maximal principal moment of inertia) form a plane. On this plane, the bone surface boundary is called 
the prime meridian, which is used as the surface cutting line to develop the bone surface. The following 
equation will make it happen: 



: x z'Vz ^ Vz 

x% = -Ax,y l z > y c z 
x% = -2Ax,yl > y c z 

4 ~ x z = 2Ax > y\ > y c z 
- x z = &x,yl > y°zi 



(6) 



where (x\, y z ) is the position of the AV at the tomographic boundary relative to the tomographic center 
of mass, (x z ,y z ) the position of tomographic center of mass relative to the inertia reference frame and i 
the sequence of AV of the boundary after being cut. 

Eq. [5] sequences the A Vs at the tomographic boundary which has been cut. The average radium of 
the tomographic boundary perpendicular to the minimal (maximal) principal moment of inertia is defined 
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as the "sea level", which is used as a datum line so that the tomographic boundary can be developed by 
the following equation: 

p(x,y) z =p(i + x c z ,hi) z , (7) 

where z shares the same definition of that in Eq. [5j h\ — r\ — r z , r\ = y 7 (x\ — x c z ) 2 + (y l z — y z ) 2 and 

r z = 2^*, (x c z ,y z ) and (x l zl y z ) have the same definitions as those in Eq. [B] 

Eqs. 0J- [7] have developed the closed surface into an open three-dimensional curved one with the 
properties of a contour. The three-dimensional curved map of the bone surface can be further developed 
into a two-dimensional plane. We can, however, translate the bone surface into a plane with the help of 
the following equation: 

P{xiy)=p(^j yp® + Z?) dt,J yp(j) + Zj) dj^j , (8) 

where p(x, y) presents the position of AV in a plane whose surface has been flattened, and 

= \J{zi- Zi-x) 2 + Zf, Zt = J \zi - Zi-i\di, = yj (zj - Zj-x) 2 + Z?, Zj = J \zj - zj^x\dj, 

z stands for the value of contour on in bone surface mapping. 

Eq. |S] suggests that the mapping of the bone surface actually serves as a simulation of the bone surface 
structure. It is a space model of an image symbol to represent the bone surface. It shares the consistency 
with the real body of the bone surface structure. 

Experiment 

Prior to our study, the Ethic Committee of Guangzhou Institute of Physical Education has proved our 
study and the participant has provided fully informed consent to participate in this study by signing a 
written consent form. From January 2008 to August 2009, we followed the track of Guangdong Provincial 
Youth Team of Male Wrestlers by using a 64 slice scanner (Brilliance 64, Philips Medical Systems). 
Excluding team members who left the team halfway and those with injuries, in January 2008 and August 
2009, we were able to scan and collect the data of a 25-year-old wrestler's sesamoid bones beneath the 
head of the first metatarsal bone of both feet. See Table 1 for information of the sesamoid bones. 

In Table [TJ the volumes of left and right foot's internal and external sesamoid bone have changed by 
0.26%, -1.98%, 0.38% and -1.32% respectively; their surface areas have changed by 0.22%, -0.41%, 
0.54% and -0.23% respectively and their density by 0.85%, 1.45%, 2.40% and 3.77% respectively. 

Let's make a mapping analysis to the right foot's external sesamoid beneath the head of the first 
metatarsal bone. First, rotate the first-time scanned sesamoid bone around axis x from to 180 degrees. 
According to Eq. [21 the variations of (I y — I z ) a , (I x — I*)fj and (I x — I v ) 1 are shown in Fig. [TJ indicating 
that within the range of rotation from to 180 degrees, an extremum exists in Eq. [2j We can get the 
result of a = 157.73 when the extremum is calculated by Eq. (3]4_. After the sesamoid bone rotates around 
axis x at 157.63 degree, it then rotates around axis y. The calculation by Eq. [3j3 generates the result 
of /? = 8.92. When the sesamoid bone rotates around axis y at 8.92 degree, the calculation of Eq. [3]C 
generates the result of 7 = 172.07. When the coordinate system of the sesamoid bone rotate around axes 
xyz at 157.63/8.92/172.07 degree respectively, the axes of of the sesamoid bone coincide with those of 
the principal moments of inertia. 

It's shown that the morphologically asymmetric and hetcrogencously distributed bone has the unique- 
ness of direction cosines of their principal moments of inertia. An arbitrary coordinate system set upon 
the bone's center of mass can make every coordinate axis coincide with the principal moment of inertia by 
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coordinate transformation, i.e. the proposed principal axis's coordinate system has its uniqueness. This 
method can be applied to both the homogeneous asymmetric geometry and the heterogeneous asymmetric 
geometry. 

The principal axes of sesamoid bone is set up by Eqs. [2] and [3J the tomography of sesamoid bone is 
reconstructed after it is rotated by Eq. [4j the boundary of the sesamoid bone is drawn by Eq. [5j When 
the prime meridian is determined by the principal axes, the sesamoid bone surface is developed by Eq. [6] 
and the mapping of the sesamoid bone surface is accomplished by Eq. Sec Fig. [2] 

Figs. [2]^ and 2B indicate that the mapping technique can better illustrate the bone surface properties. 
Fig. [5p suggests that the mapping technique can be used as a quantitative method to study the changes 
of an object's shape. Using the bone surface mapping, Eq. [8] can flatten the bone surface as a plane. See 
Fig. El 

It can be concluded that the professional training has caused adaptative changes of the sesamoid 
bone's shape. The structural changes can be depicted by the density and distribution while the shape 
and its changes can be analyzed by the bone surface mapping. 

Conclusion 

The generalized Papoulis theorem [33] elucidates that the bone surface shape keeps its geometric invari- 
ance, such as rotation, translation or dimension change [301EI] - This ensures the consistency of the CT 
scanning results of different postures of bone in vivo when its isotropy is ascertained. The uniqueness 
of the relative consistency of the inertia reference system of principal moments of inertia on direction 
cosines provide evidence to the bone surface mapping technique. The characters such as the geometric 
invariance and the uniqueness of the coordinate system of the principal moments of inertia have enabled 
the bone surface mapping technique to depict the bone's external morphological characters. This can 
advance the research of the morphological mechanisms. The experiment of the bone in vivo signifies that 
the bone mapping technique adds another research method and supplements the analytical method of 
the bone's three-dimensional imaging technique. 

We can understand the world through a map. When the bone surface mapping technique reveals 
its surface information through a "map", the activities of our life evolve continuously on this map. We 
anticipate that this mapping technique will be widely used in related disciplines. 
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Tabic 1: Sesamoid bone volume, surface area and density of two measurements 
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28.83 


210.73 


1.92 
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- internal 


205.87 


181.61 


1.89 


203.13 


181.20 


1.97 




Angular displacement 



Figure 1: The variations of the axis moment of inertia of the first measurement of the wrestler's lower 
external sesamoid bone of the first metatarsal bone accompanying the changes of right foot. 
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Figure 2: Surface mapping of bone. Fig.[2j\ Surface mapping of the first measurement of the wrestler right 
foot's external sesamoid bone beneath the head of the first metatarsal bone. Fig. [2)3 Surface mapping of 
the second measurement of the wrestler right foot's external sesamoid bone beneath the head of the first 
metatarsal bone. Fig. [2p Variations of the surface mapping of the wrestler right foot's external sesamoid 
bone beneath the head of the first metatarsal bone after 18 months (the second measurement). Figs.[5]4_, 
[2)3 and [2p have been rated by percentage and smoothed. A closed bone surface means that when it 
is continuous, there is no boundary. But when the surface is cut by the prime meridian, a boundary 
emerges. So when smoothing the surface, cloud computing method is adopted |28j to ensure the integrity 
of the object's shape. 




Principal axis of inertia (min) Principal axis of inertia (min) 



Figure 3: Flattened plane of bone. Fig. [3K The flattened bone surface of the first measurement of the 
wrestler's lower external sesamoid bone of the first metatarsal bone. Fig. [3)3 The flattened bone surface 
of the second measurement of the wrestler's lower external sesamoid bone of the first metatarsal bone 
after 18 months. 



